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Abstract 

Excitations of a vortex are usually considered in a linear approximation neglect- 
ing their backreaction on the vortex. In the present paper we investigate backre- 
action of Proca type excitations on a straightlinear vortex in the Abelian Higgs 
model. We propose exact Ansatz for fields of the excited vortex. From initial set 
of six nonlinear field equations we obtain (in a limit of weak excitations) two linear 
wave equations for the backreaction corrections. Their approximate solutions are 
found in the cases of plane wave and wave packet type excitations. We find that 
the excited vortex radiates vector field and that the Higgs field has a very broad 
oscillating component. 
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1 INTRODUCTION 



String-like vortices play essential role in several branches of physics, e.g., in condensed 
matter physics, in physics of hadrons, or in field-theoretical cosmology. The vortices have 
very interesting, intricate dynamics. Its understanding is crucial for extracting predictions 
from many theoretical models. 

We are interested in a particular aspect of vortex dynamics - properties of excited 
vortex. Excitations of vortices have been studied already for some time, see e.g. [1] 
for recent results and references to older papers. Our interest in this aspect of vortex 
dynamics stems from the fact that there is no large energy barrier or exclusion principle for 
production of excited vortices. Therefore, one may expect that in formation or interaction 
processes excited vortices will occur quite often, similarly as excited atoms or molecules 
are obtained in a recombination process. Such expectations are supported by computer 
simulations of the intercommutation of vortices [2], where one observes that parts of 
vortices which are close to the interaction point are locally excited. For concreteness, 
we consider the Abrikosov-Nielsen-Olesen (ANO) vortex [3] in the Abclian Higgs model 
with a particular type of axially symmetric excitation discussed already in [4], [5]. The 
excitation can be described as a bound state of longitudinally polarized vector particles 
with the vortex. In the present paper we study the excited vortex from a general field- 
theoretical viewpoint - applications to vortices in superconductors will be presented in a 
forthcoming paper [6]. 

The problem we are addressing ourselves to can be described as follows. The excitation 
is usually obtained as a first order correction (in an expansion with respect to amplitude 
of the excitation) to the unexcitcd ANO vortex fields, which are taken as the zeroth order 
fields. Up to the first order approximation the total field configuration is just a sum of 
the fields of unexcited vortex and of the excitation. However, because of nonlinearity of 
the pertinent field equations there are higher order corrections to this sum. Generally, 
they can be regarded as giving interaction of the excitation with the unexcitcd ANO 
vortex, and also selfinteraction of the excitation. Such second and higher order effects 
in the case of vortices seem to be unexplored as yet. In the present paper we calculate 
second order corrections in the perturbative expansion with respect to amplitude of the 
excitation in order to check what kind of physical effects appear. It turns out that these 
second order effects are quite interesting. We find that the excited vortex emits waves 
of the vector field, and that it has larger and oscillating width. All corrections found in 
the second order can be regarded as a backreaction of the excitation on the ANO vortex. 
Selfinteraction of the excitation will appear in the third order. 

The present paper is based on our earlier work [7] in which we have obtained ap- 
proximate analytic formulae for the ANO vortex, as well as for the excitation profile and 
frequency (i^^^o^ X^^*^, Oiij') and ujq below). We have also considered there the backre- 
action in the case of homogeneously excited vortex, which is a particularly simple case 
because of the lack of dependence on position along the vortex. Certain results obtained 
in [7] hold for excitations of finite amplitude, e.g., we have found an upper bound on 
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amplitude of the homogeneous excitation. In the present paper we try to cope with in- 
homogeneous excitations. The price for greater generahty is that we can investigate only 
excitations which have small amplitude - our approach is perturbative one. 

In order to keep technical details at a reasonable level we consider only elementary 
vortex with unit topological charge. The vortex is straight-linear, axially symmetric, and 
its core lies on the x^— axis. Another simplification we adopt in our paper is that the 

parameter k = \J'2q'^/ A is assumed to be very small, n <^ 1/2. In a static case this 
corresponds to strongly II type superconductors. Actually, in several places wc consider 
the limiting case k — > 0, but in fact we expect that our conclusion remain qualitatively 
valid for all K < 1/2. The point k — 1/2 is distinguished by a change of asymptotic form 
of the Higgs field [8] - a polynomial approximation used in [7] has not been extended as 
yet to this and higher values of k. 

The plan of our paper is as follows. In the next Section we present Ansatz and 
field equations for the axially symmetric excited vortex. In Section 3 we describe the 
perturbative expansion with respect to the amplitude of excitation. Sections 4 and 5 
are devoted to detailed analysis of the backreaction in the cases of plane wave and wave 
packet type excitations. In Section 6 we have collected several remarks. 



2 THE ANSATZ 

Our notation is fixed by writing the Lagrangian of the Abelian Higgs model in the following 
form, 

A / 2m'^\'^ 1 

L = {D,<p)*{Dy) - - - - i^M-^^'^ (1) 

where 

Signature of the space-time metric is (+, — , — , — ) and c = 1 as usual. Masses of the 
scalar and vector particles are equal to mn — V^m and = ktuh, respectively. In the 
following we put q = 1- 

To describe the excited vortex we use the most general axially symmetric extension of 
the Abrikosov - Nielsen - Olesen (ANO) Ansatz, 



/2m2 
^exp 



z(^ + i?(r,r))jF(r,r), (2) 

A = -niH , A = niH , (3) 

r r r r 

Af^ ^mnA'^iCr). (4) 
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Here we use the dimensionless variables, 

r'' — m.Hx'', r — ^J{r^y + (r^)^; 

The indices a, P take values 0, 3, while i,k = 1,2. 9 denotes the azimuthal angle in the 
(r^,r^) plane. The functions F, x, and ■& are dimensionless. For the unexcited, static 
ANO vortex 

= 0, = 0, (5) 

and F, x do not depend on 

The Ansatz (2-4) and Euler-Lagrange equations obtained from Lagrangian (1) lead 
to the following equations for F,'d,Al^ and x, 

.2 



- d^df^F + F" + -F' + 



1 



Q/2 



0, 



Fdpd^d + 2dfiFd^d - 2F'd' - F [d" + ^t?'^ + 2A^d0F + (a^?^^) F = 0, 
- dad'^Al^ + + (a^" + — ^ - K^F'^A'^ = K^F^d'^^, 



(6) 

(7) 
(8) 

(9) 
(10) 



Here ' stands for d/dr. Equations for the vector field we have spht into the P — 0,3 and 
i = 1,2 components. 

Equation (7) has the form of wave equation for i). It turns out that it may be omitted 
because it follows from Eqs. (8), (10). To check this, act with 9^ on both sides of Eq. (8) 
to obtain 

dpA^" = ~dpA^' + 2K^FdpF {d^d + Af^) + k^F"" {dpd^d + dpA^) . 

Next, eliminate dfsA^' from this formula using the relation (10) and divide the resulting 
formula by k^F. The result coincides with Eq. (7). 

Thus, in addition to Eqs. (6), (9) which are nontrivial also for the unexcited ANO 
vortex (cf. formulae (5)), we have two nonlinear wave equations more, namely the two 
components of Eq. (8). The phase field i) is determined from Eq. (10): 



ds F-'{C,s)d0A'^{C,s). 



This formula implies a gauge fixing for a residual gauge freedom present in the Ansatz 
(2-4) and consisting of gauge transformations with gauge function S-d dependent only 
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on Using this freedom we may adjust 'd so that it vanishes for r — cxo - this is 
the gauge imphcd by formula (11). Notice that in spite of U(l) gauge invariance of 
Lagrangian (1) the phase r) cannot be gauged away within the Ansatz (2-4) because 
the corresponding gauge transformation would introduce terms 5A^ ~ r'/r which are not 
compatible with formulae (3) for A^. In an alternative form of the Ansatz we could have 
put 'd — Q and introduced the terms ~ r^/r in formulae (3). 

Apart from the wave equations, we also require that the functions F, Xi obey a 
number of boundary conditions. For F and x at r = we take the standard ANO vortex 
conditions 

F(0) = 0, x(0) = 1, (12) 

which follow from the requirement that the fields and are regular at r = - we are 
not interested here in singular solutions. Boundary conditions at r — > oo are less obvious. 
First, in our perturbative approach the pair of equations (6, 9) is replaced by infinite set 
of equations, a pair in each order. It turns out that it is not always possible to impose 
the conditions that the perturbative contributions to F and % exponentially vanish for 
r — > oo. The reason is that most of the pertubatively obtained equations have radiation 
type solutions which vanish very slowly or do not vanish at all for large r. In this case 
we adopt the Helmholtz condition which states that for r — > oo only outgoing radiation 
waves are present, [9]. 

As for the functions A^ ^ we shall require that they are regular for all r > and that 
they vanish for r ^ oo. These conditions correspond to the physical picture that the 
excitation is localised on the vortex, i.e. that it can be regarded as a bound state of the 
/3 = 0, 3 components of the vector field with the vortex. 



3 THE PERTURBATIVE EXPANSION 

The set of nonlinear wave equations (6), (8), (9) is very complicated. In the following part 
of our paper we attempt to construct its perturbative solution relevant for the excited 
vortex in the case when the excitation, represented by the ^4^ fields, has a small amplitude 
which will be regarded as expansion parameter. The ensuing picture is as follows. In the 
leading order (e°) we have the unperturbed ANO vortex. In the order waves of the A^ 
field travelling along the vortex appear. They can be described as a 2-dimensional free 
Proca field living on the vortex. In the next order (e^) backreaction of the Proca wave on 
the vortex appears. It introduces two new features: perturbation of radius of the vortex 
and radiation from the vortex. 

The perturbative expansion has the form 

A^ = eA^^^ + 0{e^), (13) 

F = F^^O ^ ^2^(2) ^ ^^^A^^ ^ ^ ^ANO ^ ^2^(2) ^ (j^^A^^ ^^4^ 
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where F^^'-', x^^^ give the well-known static, unexcited ANO vortex. Approximate 
formulae for F^^'^, x^^'^ are given below, s is just an auxilliary book-keeping parameter 
— it is put to 1 at the end of calculations. The power expansion in e is what we mean by 
the expansion in amphtude of the excitation ^4^^^. The powers of e in formulae (13, 14) 
are suggested by the form of Eqs. (6), (8), (9). 

Let us first eliminate the phase In the leading order, Eq. (10) admits as the solution 

= 0, (15) 

with 

d^A^fj^ = 0. (16) 

This solution is compatible with Eq. (8) because d/sF^^^ = 0. 

Now, with (15) taken into account, Eq. (8) in the order decouples from the other 
equations, 

- + 4'^" + ^4')' - K^F^'^yA^fj^ = 0. (17) 

This equation has the form of 4-dimcnsional wave equation and as such it has plenty of 
solutions. We are interested in bound state type solutions, i.e. the ones which exponen- 
tially vanish for large r. They can be found with the help of separation of variables 
from r, 

4^^W^{0a{r), (18) 
where VF^(^) obey the Lorentz condition 

d^VV^ = 0, (19) 

and the 2-dimensional Proca equation 

d^d'^Wfs + ujlWfs^O (20) 

with the mass cuq given below. 

It follows from (17-20) that the profile function of the excitation, i.e. a{r), obeys the 
equation 

u;o'«(r) + a"(r) + ^ - «2(F^^On2^/ ^ _ q (21) 
r 

This equation has been considered in [5], [7]. Its vanishing for r — > cxo solution is quoted 
below. Without any loss of generality we may normalize a{r). The solution given below 
is normalized by the condition 

a{0) = 1. 

Proca equation (20) possesses infinitely many solutions. Among them there are trav- 
elling plane waves and wave packets - they give rise to the inhomogeneous excitations 
considered in the present paper. The "breathing" homogeneous vortex considered in [7] 
is obtained for W'^ = 0, = W^{x^), which is a non-propagating, periodic in time and 
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x^-indepcndcnt Proca field standing on the vortex. Condition (19) implies that there 
exists a potential ip for the field, 

where e^Q, = — e^/j and eo3 = +1. Proca equation (20) will be satisfied if the potential 
obeys Klein-Gordon equation 

+ coli/j = 0. (22) 
In the following Sections we shall consider the plane wave solution of equation (22), namely 

,p = Nsm{E{q)e-qe + So), (23) 



where E{q) = \JuJq + q^, N is a constant, and Sq is a constant phase which we shall put 
to zero. 

Another solution of Eq. (22), also considered in the following Sections, is an approxi- 
mate solution giving a standing wave packet, 

i; ^ Nexp{-ej^A') sin(e%), (24) 

where A and N are constants. Average momentum of this wave packet is equal to zero. 
Wc also assume that ujqA ^ 1. Then, the wave packet has momentum cutoff at small 
momentum ~ A^^ <^ 1. In this case we may neglect (for a finite time < < A^) 
dispersion of the wave packet, which is of course present in exact wave packet solutions 
of Eq. (22). Such exact wave packet solutions are easy to obtain, nevertheless we prefer 
the approximate solution because it has significantly simpler form. 

Now, let us turn to the remaining two equations (6) and (9). In the order they are 
solved by the unexcited vortex fields F^^^, x^^'^- In the order s"^ they give the following 
equations for x^'^\ ^^'^^ '■ 



r 2 



1 - 3(F^^O)2 _ ^(^ANO)2- 

2 



(25) 



r 



2/ 



_ df,d^x^^^ + X^'^" - ^X^'^' - «2(F^''°)'x^'^ = 2/^2f^nOx^''°^^'^- (26) 

The boundary conditions for F^'^\x^'^^ at r = follow from conditions (12) and from the 
fact that i?ANO -j^ANO already obey these conditions. Therefore, 

F(2)(r = 0,r) = 0, x^'\r = 0,C) = 0. (27) 
We shall see that we may also require that 

^^^H^D^O for r^oo. (28) 
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As for x^'^\ we find that Eq. (26) has radiation type solutions, so we adopt the outgoing 
radiation condition. 

The last term on the r.h.s of Eq. (25) is a source term for F^'^\ Due to its presence 
^(2) T^(2) jjq-j- vanish. Hence, the backreaction of the excitation on the vortex is always 
present (one can easily prove that WpW^ 7^ if Wp ^ 0) . 

Prom a mathematical point of view, Eqs. (25), (26) form a set of 4-dimensional, linear 
wave equations with the coefficients given by the ANO functions ^t^ano^t^ano gofiitions 
of such equations can be investigated with a help of numerical methods. In the present 
paper we take another route. We will present approximate analytical solutions which we 
can obtain in the limiting case k — > 0, i.e. when < k <S 1/2 (in practice we mean k not 
greater than 0.1). This limit is distinguished by two significant simplifications. 

First, because /t < 1/2 we may use relatively simple approximate formulae for the 
functions p^^'^ ^ X^^*^, OL{r^ and for the frequency cjq obtained in [7]. In the limiting case 
K <^l/2 we obtain from formulae given in [7] that for r < ro = 4/^3 



while for r > Tq 



where 



and 



X*"" ^ 1, (29) 

a{r) ~ 1, (31) 



X^""^ ^ KrXi(Kr), (32) 

F^^o ~ 1, (33) 

a{r) ~ ciKoikor), (34) 

23 o . . 

Ci ~ , (35) 

A;o~0.85 exp{-^^}. (36) 



The frequency cuq is given by the formula 

uo = \Jk'^ - kl- 

Kq[z) and Ki[z) are the modified Hankel functions, [10]. The corresponding functions in 
the two regions smoothly match each other at r = Tq. On the r.h.s. of formulae (29-36) 
we have neglected terms of the order In k or smaller. More accurate formulae for the 
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ANO vortex and the excitation can be found in [7]. Notice that ko is exceedingly small, 
and therefore the profile function a{r) is very broad, much broader than x^'^^^ and F^^^ 
which have the width ~ I/k, ~ Tq, respectively. 

The second simplification occuring for k <^ 1/2 is that the first term on the r.h.s. 
of Eq. (25) can be neglected. To see this, first notice that this term is finite at r = 
because F^^*^ and x^^"* vanish at that point. Now, the r.h.s. of Eq. (26) is of the order k^, 
therefore x^"^^ determined from that equation is expected to vanish for k ^ 0. Therefore, 
the discussed term for the very small k will be negligibly small in comparison with the 
other term on the r.h.s. of Eq. (25) which stays finite in the limit k — > 0. Thus, instead 
of Eq. (25) we may consider the equation 



Due to this latter simplification the set of intercoupled equations (25), (26) is replaced by 
equations (37) and (26) which can be solved one after the other. F^^^ determined from 
Eq. (37) gives the source term in Eq. (26). In the following two Sections we shall present 
solutions of these equations in the particular cases of the excitation field given by the 
plane wave (23) and the wave packet (24). 



4 BACKREACTION ON THE HIGGS FIELD 



In this Section we shall analyse the backreaction of the excitations on the Higgs field F. 
Our starting point is the equation (37). 

It is convenient to pass to Fourier transforms with respect to 



+ F(2)" + 1f(2)' + 1f(2) [l - 3(F^^O)2 _ 2 (^ANO)2 



(37) 




(38) 




(39) 



and 




(40) 



The Fourier transform of Eq.(37) has the following form 



f" + -f -\oo'-k'-U |(F^^°(r))^ + ^(X^^°(r))^l / 
r L 2 2 r 



(41) 



F^''°{r)a''{r)w{k,u;). 



For the plane wave (23) with 6o = 



8 



and 

w{k, uj) 



6{ij)6{k) + -6{u - 2E{q))5{k - 2q) + -5{uj + 2E{q))S{k + 2q) 
2 2 



In the case of the wave packet (24) 

W^Wa = w.iis) + ^2(^3) cos(2^oe°), 

where 



el 



2A2^' 



^2(e3) = -7,A^^(^+c.o^)exp(-;^). 



2A2' 



Formula (40) gives in this case 
w{k, u) — V27r 



5{u})wi{k) + ^S{u} - 2u}o)w2{k) + ^S{u + 2uo)w2{k) 



(42) 
(43) 

(44) 
(45) 

(46) 



where 2 denote Fourier transforms of ^1^2(^3)- We shall not need their explicit form. 
It is clear from Eq.(41) that / has the form analoguous to (42), (46), respectively. 



fo{r)6{uj)6{k) + -U{r)S{uj - 2E{q))5{k - 2q) 



+-f4r)5iu; + 2E{q))5{k + 2q) 



(47) 



or 



fo{r)S{u;)wi{k) + lf+{r)S{uj - 2uJo)w2{k) + \f-{r)8{uj + 2uJo)w2{k) 



■ (48) 



The negative and positive frequency components are related by complex conjugation, 

/- = /; 

and /o is real valued. 

Next, we observe that the expression in square bracket on the l.h.s. of Eq. (41) can 
be simplified a little bit. Namely, the term cj^ — k'^ is equal to either or AuJ^ in both the 
plane wave and wave packet cases, and for small k it is negligibly small in comparison with 
the other terms in the square bracket - a plot of the function -1/2 + 2,/2{F^^^{r)f + 
(x^^*^(r))^/r2 shows that its minimal value is approximately equal to 0.945. Therefore, 
instead of Eq. (41) we may consider the following universal equation, common for all 
frequency components and for the plane wave and wave packet cases. 



f" + -f' - 
r 



fu 



nANO, 



r)Q;^(r), 



(49) 



9 



where 

/« = /o ~ f±. (50) 

In order to see the backreaction on the Higgs field we have to solve Eq. (49) for /„. 
To this end we shall use the polynomial approximation which has proved useful also in 
our earlier work [7] . The idea is to find approximate solutions of Eq. (49) separately in 
the regions r < Tq and r > Tq, and to glue them at r = Tq. 

In the "outer" region, i.e. for r > ro — 4/-\/3, we use formulae (32-36). Thus, 
F^^'^{r) ~ 1, and the functions (-i^ano-)2^^2 a{rY are smooth and very slowly chang- 
ing with r. Their derivatives with respect to r are of the order k and ko, respectively, hence 
they are close to zero. Therefore, the following function is a reasonable approximation to 
the exact solution of Eq. (49) in the outer region 

/„ = /„(r)^-a(r)^(^l+ ^^ j + foMr), (51) 

where the last term on the r.h.s. is a general, vanishing for r — > oo, solution of the 
homogeneous counterpart of Eq. (49). /o is a constant to be determined from matching 

conditions at r = Tq. 

On the other hand, for r < tq the ANO functions are given by polynomials (29-31), 
and we seek a polynomial approximation also for /„, 

f^^l^ = f,r + fsr' + fr,r'. (52) 

We take the fifth order polynomial for /„ because formula (30) gives F-^^'^ up to the 
term. Inclusion of the next term (~ r^) would make sense if we knew F^^'~' up to ~ 
term, as it is clear from comparison of the both sides of Eq. (49). Equation (49) gives 
the following relations 

■'^ 16ro le-' ' 
19 1 

^ ~4S^ ^ e^f-^' ~ 48-^'' 

where /i remains undetermined. 

At r = To we impose matching conditions. They have the following form 

/„(ro) = /.(ro), £(ro) = /^(ro). 

These conditions give a set of linear, algebraic equations for the constants /o, /i - the 
solution is 

/of«5.19, fi^-OM. (53) 

From formulae (38), (47), (48), (50) we obtain the following final expression for the 
backreaction on the scalar field 

^ -WW^iOfuir), (54) 
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where /^(r) is given by (51-53). Formula (54) covers the plane wave and wave packet 
cases. We see that F^^^ introduces modulated along the vortex and oscillating in time 
very broad component in the scalar field of the vortex. Its range is ~ (korriH)'^- 



5 BACKREACTION ON THE VECTOR FIELD 



Now let us investigate the backreaction on the vector field. Fourier transform of equation 
(26) has the following form 



1./ 



h"--h'+ uj^-e-^{F^''''{r)f /i = 2/^2F^^"(r)x^^"(r)/(r,A;,u;). (55) 
We shall consider the plane wave and the wave packet cases separately. 



5.1 The plane wave case 

In analogy to formula (47) we write 



h = 2k\ujIN^ 



ho{r)5{uj)5{k) + -h+{r)6{u - 2E{q))6{k - 2q) (56) 
+hi^{r)5{u + 2E{q))5{k + 2q) 

Equation (55) is equivalent to the following set of equations (one equation for each 
frequency component), 



"•a 

r 



K - ^\F^''^{r) f K = F"^^"(r)x"""(r)/„(r), (57) 



'ANO/„N,,ANO/ 



where the index a takes values 0,+,—, and VLq = 0, Q± = ±2uJo. Very important 
difference with the case of Higgs backreaction / is that in these equations the 4uJq term 
is of the same order of magnitude as the k,'^{F^^^)^ term (recall that Uq ^ /t^). 

There are two essentially different cases: fl^ = and flj_ = AuJq. In the first case, in 
the outer region (r > Tq) Eq. (57) has the form 

K-h',-K%^x''''''{r)Ur), (58) 

where x^^° is given by formula (32). Applying standard methods, see e.g. [11], we find 
the following exact, vanishing at r ^ oo, solution 

ho{r) — rgi{r)Ii{Kr) + rg2{r)Ki{Kr) + boKrKi{Kr), (59) 

where 

rfr'x^''°(r')i^i(«:r')/„(r'), 92{r) = - dr' x''''''{r')h{Kr')fu{r'). (60) 

'ro 
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Here Ji is the modified Bessel function [10]. The constant bo will be determined from 
matching conditions at r = Tq. We see that for very large r, i.e. for r ^ /cq ^, 

ho{T^) ~ -\/r exp(— Acr). 

This component of the backreaction on the vector field is relatively uninteresting. It gives 
a small, time- and ^^-independent, correction to the t^^no^^^ ^j^g same range as 

^ANO ^^jjg penetration depth ~ «;~^m^^). 

Much more interesting are the other components, i.e. the ones with Q,± — ±2a;o. In 
this case, Eq. (57) reduces in the outer region to 

h'^ - h'^ + (4^o' - ^'yh± = X^^°(r)A(r), (61) 

where again X'^^^ is given by formula (32). Applying the standard methods [11], we 
obtain the solution 



= rhl{r)H'^\^J AujI - k"^ r) + rh2{r)H'^\^J AujI - k"^ r) 

(62) 

+6±^4cu2 - rH^^'^\^jAujl- r), 
where h\^\h^^^ = (i/f^) * are Hankel functions, and 

h,{T) = I dr' Hf\^Au^l-K^ r')x^^°(r')/«(r'), (63) 

In order to keep x real-valued we have to assume that 6+ = (^-)*- Iii the last term on 
the r.h.s. of formula (62) one should take Hi^^{h'i'^) in the positive (negative) frequency 
component in order to satisfy the outgoing radiation condition. 

The functions hi^2 exponentially vanish for large r. On the other hand, for large z [10] 



V TTZ 



exp(i(^-37r/4)), (64) 



hence the last term on the r.h.s. of formula (62) grows like r^^^. It gives the radiation of 
the vector field from the vortex. 

In the small r region (r < ro), for the all frequencies, we shall be satisfied with 
approximate polynomial solutions in the form 

^-4^-^^+^^^+J^(^3-^).^ (65) 

h_ = {h^y. 
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They have been obtained by inserting the formulae (29, 30, 52) in Eq. (57), and by 
neglecting all terms ~ on the l.h.s. of it (let us remind that we consider the case 
K <^l/2). Again, we have taken a polynomial of the maximal order compatible with the 
fact that for small r the r.h.s. of Eq.(57) is known up to terms of the order r^. 

The constants bo,b±,do, d± are determined from matching conditions at r = tq: 

ho,±{ro) = ho,±{ro), /i|),±(ro) = /io,±(ro). 

It turns out that the first two terms on the r.h.s. of formula (62) computed for r = tq are 
of the order k, ln(K) and therefore they can be neglected for small k. The argument of the 
Hankel function Hj^^ is very small for r = Tq, so we may use the following approximate 
formula [10] 

Hi\z) ~ -iOM/z. 

Similarly, 

KToKi^nro) ^ 1. 

The resulting from the matching conditions sets of linear algebraic equations, simplified 
as described above, give 

13 5 2 

do ~ Kgi{ro) - -/iro - -, 6o ~ "g/i^o - g, (66) 

d± ^ -Ifir^' - 2ro ^ (67) 

and 

b± ^ ±1.56i6o T -inK'\iujl - k^)-^^'^C{k), (68) 
where C{k) is a constant given by the following formula 

C(«)= r ds Ji(^3 - 4A;g«-2 s)sK^{s)ju{slK). 

In formula (68) only the second term on the r.h.s. is important because is very large. 
The constants tq and /i we know from Sections 3 and 4. This completes the determination 
of the backreaction on the vector field. 

The most interesting implication of these formulae is that the backreaction includes 
radiation from the excited vortex. The b± components of h± inserted in formulae (56), 
(39) give the following contribution to x^^^ for r » Tq 

Sx^^^ fti -KwliAul - n^)-^/^N^^/^C{K)V?cos{-2E{q)^'' + 2q^^ + Aul - r - 7r/4), 

(69) 

where E{(i) — + q^. This radiation wave has a conical shape - at fixed time the 
surface of constant phase, determined by the following condition 

2qi^ + sj^ul - r^T:/A + 2E{q)^'^ + const, 
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is a cone with its top point lying on the axis. For q = the cone degenerates to a 
cyhnder of constant r. Opening angle 70 of the cone, given by the formula 

2q 



tan 7o — 



increases with the wave number q of the Proca wave on the vortex. The wave four-vector 
corresponding to the radiation wave (69) has the following components 

{kn^{2E{q),k), 

where 

k — \J4wq — (r^/r, r'^/r, tan 70) . 

Its 4-dimcnsional invariant length is equal to k^, in accordance with the fact that the 
vector particle has the mass equal to k in the uih units. 

The energy flux related to that conical wave is given by the Poynting vector 

S = 5Ex 6B, 

where 5E, 5B are the contributions to electric and magnetic fields corresponding to 5%*^^^. 
Simple calculation gives 

S = 7rC(«)2«2,^^(4^2 _ ^Y'^^ N^rn%E{q)^ sin'' {k^x^mH)k. (70) 
5.2 The wave packet case 

Calculations of backreaction of the Proca wave packet on the vector field x ^-^e carried 
out in the same steps as in the plane wave case. Instead of formula (56) we have now 



/i = -2(27r)VV 



wi{k)5{uj)ho{r) + ^W2{k)S{uj - 2ujo)h+{r) (71) 



Equations for the frequency components /iq, ^± coincide with Eq. (57), where is 
replaced by — fc^ or Alo^ — k'. However, taking into account the assumption that the wave 
packet (24) has momentum cutoff at /c ~ where oo^h. ^ 1, we see that for our wave 
packet k'^ <^ oj^ so that we may neglect the —k'^ term altogether. Hence, we again obtain 
Eq.(57) with Vt^ = 0, Vt\ = Auq as in the plane wave case. Therefore, the discussion 
following Eq. (57) apphes also to the present case, including results (66-68) for constants 
d±, c±. We again find the radiation. Using formula (62) and asymptotic formula (64) for 
the Hankel function we obtain the following form of the radiation part of the backreaction 
in the standing wave packet case 

Sx^'^ ^ (27r)^/2«(4u;2 - «')-'/'C(«)«;2(e3)v^cos(-2a;oe' + 2?^' + ^^^o - r - 7r/4), 

(72) 
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It is a cylindrical wave, modulated along the vortex by the ^2(^3) function. The corre- 
sponding wave vector has the components 



1 1 

— T 




Its 4-dimensional invariant length is equal to k^. 

The fact that the wave front has the cylindrical shape is in accordance with the conical 
character of the emitted wave in the previous example. The standing wave packet can be 
regarded as containing +q and —q components with an equal weight - such a superposition 
of conical waves (69) gives the cylindrical wave. 

Proca wave packets moving along the vortex can be obtained by Lorentz boosts. 



1. Let us summarize the main points of our paper. In order to investigate the excited 
vortex fields, introduced by the Ansatz (2-4), we have used the perturbative expansion in 
the amplitude of the excitation. The backreaction is given as the corrections F'^'^\ x^"^^ to 
the unexcited vortex fields F^^'^^ X^^^- We have found that F^"^^ exponentially vanishes, 
but the corresponding characteristic length is large (~ kQ^m'jl), much larger than the 
coherence length (~ rn'jl) and the penetration depth (~ K''^m'jl). For the very small 
values of k assumed in our paper, there is no radiation of the scalar field from the excited 
vortex. The vortex fields F, x couple to WaW" - this composite field has an effective 
mass equal to Auo^. For /t ^ 1/2 this characteristc frequency turns out to be too small 
to induce radiation of the Higgs field F, but it is large enough for radiation of the vector 
field X- Nevertheless, radiation of the Higgs field might appear in higher orders in the 
perturbative expansion. The frequency of the radiative component of 5x^'^\ c.f. formulae 
(69), (72), is twice of that of the excitation field Wa{^). 

2. Perhaps the main question about our approach is whether the approximations used 
are correct. Rigorous convergence proofs arc not known neither for the polynomial ap- 
proximation nor for the expansion in the amplitude of the excitation. Confrontation of 
the polynomial approximation for the unexcited ANO vortex with a purely numerical 
solution shows that it works rather well [12]. One reason is that the polynomials are 
used in a finite interval of r (0 < r < ro), where the vortex fields are expected to have a 
simple behaviour (e.g., one docs not expect that they oscillate like rsin(l/r) for r — >■ 0). 
Therefore, low order polynomials can give quite good approximation. 

Another reason is that the scalar and vector fields in the Abelian Higgs model are 
massive ones, and therefore the vortex fields approach the corresponding vacuum values 
very rapidly (exponentially). Therefore, deviations from the vacuum values for large r 
can be regarded as small corrections which one may reliably determine from linearised 
equations. 
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REMARKS 
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As for the expansion in powers of the amphtude, certainly encouraging is the fact that 
the correction F^'^^ has turned out to be small for all values of r. So is x*^^V^ (which is 
the function actually present in our formulae). Moreover, the exact equations (6,8,9,10) 
depend on the excitation fields A^, i? in a smooth manner, and the corresponding terms 
are just additional contributions to the other nonlinear terms which arc already present 
in the case of the imexcitcd ANO vortex. The presence of the excitation fields in 
Eq. (6) does not seem to introduce any singularity in that equation. All that suggests 
that also the solution describing the excited vortex depends smoothly on the amplitude 
of the excitation field if the amplitude is small. 

3. The presence of radiation means that the excitation is not stable. Nevertheless, 
lifetime of the small amplitude excitations is large because energy of the excitation is 
proportional to A^^ (for concretencss we arc referring to the plane wave excitation) while 
the energy fiux is of the order A^'^. Our formulae also show that the energy fiux increases 
when the amplitude grows, however our approximations should not to be trusted when 
the amplitude becomes large. 

4. Finally, let us mention two directions for rather interesting extensions of the present 
work. According to our results, the excited vortex is no longer localised within the pene- 
tration depth range ((«;mjj)~^). This likely infiuences its interaction with other vortices. 
The new interaction, unknown as yet, may turn out to have quite interesting properties. 

Also, it would be interesting to apply the polynomial approximation in order to cal- 
culate radiation of Goldstone particles from global vortices. Such an analytical approach 
would nicely complement existing numerical investigations as well as analytical results 
based on Lund-Regge string model for the global vortices presented in, e.g., [13], where 
one can also find references to older papers on radiation from the global vortices. 
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